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摘 要 : 针对 经 验 模 态 分 解 (empirical mode decomposition，EMD) 过 程 中 存在 的 包 络 拟 合 问题 ， 提 出 了 一 种 消减 欠 冲 
现象 的 改进 算法 。 该 算法 通过 引入 “ 伪 极 值 点 ”， 增 加 了 极 值 点 的 数目 ， 构 成 了 新 的 极 值 序列 ， 然 后 利用 新 的 极 值 序列 
插值 拟 合 得 到 新 的 包 络 线 。 最 后 通过 仿真 实验 对 比 本 文 算 法 和 经 典 拟 合 算法 包 络 拟 合 产生 的 欠 冲 点 数目 。 实 验 结果 显 
示 ， 与 经 典 拟 合算 法 相 比 ， 改 进 的 算法 产生 的 欠 冲 点 数目 减少 了 大 约 77.5%。 实 验 结 果 表 明 ， 此 算法 可 以 有 效 地 消减 
欠 冲 点 的 数目 ， 拟 合 出 的 包 络 线 更 加 贴近 原始 信号 ， 拥 有 更 好 的 平滑 
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Abstract: Due to the envelope fitting problem exists in the process of empirical mode decomposition, this paper proposed an 
improved algorithm which could eliminate the undershoot phenomenon exactly. By introducing "pseudo-extreme points", this 
algorithm increased the number of extreme points and formed new extreme value sequence. Then it got new envelope by using 
the new extreme value sequence interpolation. The envelope fitted by this method was closer to the original signal and has better 
smoothness. Finally, A contrast result between cubic spline interpolation and this algorithm showed that the number of 
undershoots decreased by about 77.5%, which proved that this algorithm could effectively reduce the number of undershoot 
points. In addition, fitting the envelope could tightly wrap the original signal and have a better envelope. 
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0 引言 这 些 方法 在 一 定 程 度 上 改善 了 EMD 分 解 过 程 中 包 络 线 的 拟 合 
村 性 ， 但 在 处 理 非 线 性 、 非 平稳 信号 时 ， 包 络 线 的 拟 合 特性 没 
让 电信 号 是 一 种 非 平稳 、 非 线性 的 随机 信号 ， 单 独 的 考虑 有 了 明显 改善 ， 依 然 出 现 了 许多 欠 冲 现象 。 
时 域 特 征 或 频 域 特征 都 是 不 全 面 的 ， 所 以 越 来 越 多 的 研究 致力 本 文 针 对 EMD 分 解 非 线 性 、 非 平稳 信号 的 包 络 拟 合 问题 ， 
于 脑 电信 号 的 时 频 特 征 的 研究 。EMD 是 Huang 等 提出 的 一 种 “以 脑 电信 号 为 例 ， 着 手 ， 提 出 了 一 种 消减 
自 适应 信号 的 时 频 分 析 方 法 ， 它 能 把 复杂 的 数据 分 解 成 若干 个 ” 欠 冲 现象 的 改进 方法 。 该 方法 通过 引入 “ 伪 极 大 值 点 ”和 “ 伪 
有 限 的 本 征 模 态 函 数 (IMF), 分 解 的 IMF 能 够 体现 信号 的 非 线 极 小 值 点 ”， 以 此 扩大 极 值 点 序列 , 增加 极 值 点 数目 , 来 消减 欠 
性 、 非 平稳 性 特征 ("1。 因 为 信号 的 极 值 点 序列 的 包 络 线 直接 影 。” 冲 现象 。 同时， 由 于 极 值 点 数目 的 增加 ， 也 提高 了 拟 合 曲 线 的 
响 着 EMD 分 解 的 收敛 速度 和 准确 度 ， 所 以 信号 包 络 线 的 拟 合 ”平滑 性 。 
外 在 EMD 分 解 过 程 中 至 关 重 要 。 1 ”经 典 EMD 算法 的 原理 
目前 ,研究 人 员 已 经 提出 了 多 种 包 络 拟 合 改 进 的 方法 :B 样 
条 插值 中 、 分 段 贤 函数 插值 外 、 分 段 抛物 线 插值 中 、 分 段 三 次 信号 进行 EMD 分 解 后 ， 将 会 获得 一 系列 的 频率 范围 不 同 
Hermite 插值 站 、 最 小 长 度 约束 上 中、 二 阶 包 络 逆 EMD [0 等 等 。 的 IMF 分 量 ， 每 一 个 IMF 分 量 必须 要 满足 两 个 条 件 帆 ，a) 在 获 
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取 IMF 信号 的 过 程 中 , 信号 的 局 部 极 值 点 的 个 数 和 过 零点 的 数 
目 相等 或 最 多 相差 1; a 极 大 值 构 成 的 上 包 
络 线 和 由 极 小 值 构 成 的 下 包 络 线 的 均值 为 0。 

EMD 分 解 的 具体 过 程 11- 引 如 下 : a) 求解 原始 信号 x(?) 的 
所 有 的 局 部 极 大 值 点 和 极 小 值 点 ，b)〉 对 所 求 的 极 大 值 点 和 极 
小 值 点 分 别 进行 三 次 样 条 插值 ， 求 出 上 包 络 线 xu(D 和 下 包 络 线 


zt(0e) 计 算 上 下 包 络 线 的 均值 册 (D= 玖 人 二 人 .4) 将 原始 信 


号 x( 四 减 去 包 络 均值 mi 得 到 一 个 去 掉 低 频 的 信号 

及 (四 =x( 四 一 m (四 ;e》 判 断 所 求 出 的 加 (DD) 是 否 满足 IMF 定义 的 
两 个 条 件 ,如 果 满 足 ， 户 (0 就 是 筛选 出 的 第 一 个 IMF 分 量 ， 如 果 
不 满足 ， 则 万 (D 代 替 原 始 信号 x(?)， 重 复 上 述 过 程 有 次， 直到 
生 (1) 满 足 IMF 分 量 的 条 件 44(D) = (2?)= Ry (DD) 一 mx (DD) ;f) 
用 原始 信号 x() 减 去 c1(?) 得 到 一 个 去 掉 高 频 成 分 的 信号 m0)( 残 


王 成 龙 ， 


高 了 拟 合 曲线 的 平滑 性 。 有 具体 过 程 如 下 : 
a) 提取 原始 信号 x(0 的 所 有 的 局 部 极 大 值 点 和 极 小 值 点 。 
b) 计算 相 邻 数据 点 之 间 的 斜率 。 并 计算 原始 信号 数据 点 之 
间 的 斜率 之 比 。 


d(x(D) 
dt 


c) 由 于 原始 信号 数据 点 为 极 值 点 时 ， 和 斜率 之 比 小 于 零 。 并 
且 前 后 斜率 一 致 时 ， 和 斜率 之 比 为 1， 所 以 将 计算 后 的 斜率 之 比 
减 去 1， 得 到 d' (反映 曲线 的 弯曲 程度 )。d' 的 值 越 大 ， 说 明 原 
始 信号 数据 点 前 后 越 弯 曲 。 如 果 d 大 于 0 并且 d' 大 于 0.5, 则 判 
定 该 数据 点 为 “ 伪 极 值 点 ” 记 为 D(m)。 
D(m)=d>0&d'’>0.5 


差 ) (1)=Xx(1) 一 G();g) 将 ri(D) 作 为 新 的 原始 信号 


， 重 复 上 述 


d) 对 所 求 出 的 “ 伪 极 值 点 ”进一步 划分 为 “ 伪 极 大 值 点 ” 


过 程 ,直到 循环 n 次 得 到 的 残 差 为 一 个 六 


和 调 函 数 或 是 常量 为 止 ， 


和 “ 伪 极 小 值 点 ”如果 原 始 信 


号 数据 点 和 原始 信号 数 ] 


时 点 前 后 


因此 原始 信号 经 过 EMD 分 解 后 得 到 x(1) = AOAGE 


原始 的 三 次 样 条 插值 


一 一 一 原始 信号 
上 包 络 线 


欠 冲 区 间 __、 


下 4 1 L 1 1 L 1 L 
50 55 60 65 70 75 80 85 390 95 100 


图 1 包 络 线 拟 合 的 欠 冲 现象 


从 上 述 过 程 可 知 ， 包 络 线 拟 合 贯穿 着 整个 EMD 分 解 的 过 
程 。 但是， 经 典 的 EMD 使 用 的 三 次 样 条 插值 拟 合 出 来 的 包 络 
线 ， 不 能 有 效 的 包 里 原始 信号 ， 导 致 出 现 欠 冲 和 过 冲 现象 ， 进 
而 会 产生 错误 的 极 值 点 ， 影 响 原 有 极 值 点 的 位 置 ， 造 成 较 大 的 
EMD 分 解 误差 ， 使 得 EMD 分 解 失去 意义 9。 


2 ， 包 络 线 拟 合 的 改进 


为 了 更 直观 地 展现 经 典 EMD 算法 包 络 拟 合 过 程 中 ， 出 现 
的 欠 冲 现象 ， 以 一 段 来 自 Koelstra 等 人 07 提 取 的 脑 电 情感 分 析 
数据 库 中 的 脑 电信 号 为 例 。 经 典 拟 合算 法 进行 包 络 拟 合 过 程 中 
出 现 的 欠 冲 现象 如 图 1 所 示 。 产 生 这 一 现象 的 原因 0 有 两 个 : 
a) 进 行 包 络 拟 合 时 只 考虑 了 极 值 点 的 位 置 ， 没 有 对 极 值 点 之 间 
的 斜率 进行 约束 ; b) 找 出 的 极 值 点 的 位 置 并 非 是 真正 的 极 值 点 
的 位 置 。 针 对 上 述 现象 ， 本 文 针 对 EMD 分 解 非 线性 、 非 平稳 
言 号 的 包 络 拟 合 问 题 ， 以 脑 电 信号 为 例 ， 提 出 了 一 种 消减 欠 冲 
现象 的 改进 方法 。 该 方法 通过 扩大 极 值 点 的 范围 ， 增 大 极 值 点 
数目 ， 来 消减 欠 冲 现象 。 同 时 ， 由 于 极 值 点 数目 的 增加 ， 也 提 


构成 的 曲线 是 凸 的 , 则 该 数据 点 为 “ 伪 极 大 值 点 ” 如果 原 始 信 
号 数据 点 和 原始 数据 信号 数据 点 前 后 构成 的 曲线 是 止 的 ， 则 该 
数据 点 为 “ 伪 极 小 值 点 ”。 


x(D(m-1))+x(D(m+1)) 
2 


max =x(D(m)) > 


x(D(m-1))+x(D(m+1)) 
2 
e) 将 求 出 的 “ 伪 极 大 值 点 ”和 “ 伪 极 小 值 点 ”与 极 大 值 点 
和 极 小 值 点 分 别 结合 ， 构 成 新 的 极 值 序列 。 
于 扩大 了 极 值 序列 的 范围 ， 采 用 改进 的 算法 对 信号 x(?) 
进行 包 络 拟 合 ， 拟 合 的 曲线 更 加 贴近 原始 信号 。 结 果 如 图 2 所 


min =x(D(m))< 


| 


10 玉 Tt - - — - ~ - + 

原始 信号 

a 三 次 样 条 插值 
一 一 一 改进 包 络 线 


J 1 1 hn 1 1 1 h h J 
50 55 60 65 ?0 75 80 85 90 95 100 


图 2 改进 前 后 包 络 线 拟 合 的 对 比 图 


3 ”改进 算法 仿真 结果 分 析 


为 了 验证 本 文 提出 的 改进 的 包 络 拟 合算 法 的 有 效 性 ， 以 8 
064 个 采样 点 的 脑 电信 号 为 例 。 将 本 文 包 络 拟 合 的 改进 算法 与 
经 典 包 络 拟 合算 法 、 分 段 抛物 线 插值 拟 合算 法 和 三 次 Hermite 
插值 拟 合算 法 ， 进 行 了 对 比分 析 。 具 体 结果 如 图 3 所 示 。 


201804.02156v1 


chinaXiv 


录用 稿 


400F T T T T T T 


一 经 奥 拟 合 一 分 段 扫 物 线 拟 合 一 三 次 Hermite 拟 合 一 本 文 算法 拟 合 


0 1 | 1 | 1 
0 1000 2000 3000 4000 5000 6000 


其 中 ， 图 中 横 坐 标 表示 采样 点 数目 ， 纵 坐标 表示 欠 冲 点 数目 。 


1 
?000 8000 


图 3 相同 采样 点 下 几 种 包 络 拟 合 算法 欠 冲 点 数目 的 变化 


表 4 采样 点 为 8000 时 不 同情 绪 标 签 的 欠 冲 点 数 
算法 love exciting depressing hate 
经 典 拟 合 353 355 355 374 
分 段 抛物 线 206 212 194 197 
三 次 Hermite 181 192 173 174 
本 文 算法 105 59 65 90 


表 中 可 以 看 出 ， 本 文 算法 显著 的 减少 了 欠 冲 点 数目 。 当 


采样 点 数 分 别 为 2000、4000、6000 和 8000 时 ，4 种 不 同情 绪 


EMD 分 解 过 程 中 ,使 用 经 


包 络 把 
中 蓝 色 线 条 ) 进行 拟 合 ， 产 生 的 欠 冲 点 数目 ， 随 着 采样 点 数目 


合算 法 (图 


标签 


使 ) 


j 本 文 算法 对 脑 电信 号 拟 合 


产生 的 欠 冲 点 数目 之 和 ， 


与 分 段 抛 物 线 插值 拟 合算 法 产生 的 欠 冲 点 数目 相 比 ， 分 别 少 了 


96 个、 


点 数目 近似 呈 直 线 上 升 趋势 。 


与 经 


包 络 拟 合 算 


的 增加 ， 欠 六 


Heimite 插值 拟 合算 法 (图 中 绿 


比 ， 分 段 抛 物 线 插值 拟 合算 法 (图 中 红色 线条 ) 和 三 次 
色 线 条 )， 在 一 定 程度 上 减少 了 


欠 冲 点 数目 。 而 使 ) 


行 拟 合 ， 随 着 采样 点 数 
为 8000 时 ， 经 典 包 络 拟 合 


包 络 拟 合算 法 产生 的 采样 点 数目 的 4.3 倍 。 
为 了 进一步 验证 本 文 算法 的 鲁 棒 性 ， 采 用 脑 电 情感 数据 库 


改进 的 包 络 拟 合算 法 《图 中 黑色 线条 ) 进 
的 增加 ， 上 升 趋势 缓慢 。 
算法 产生 的 采样 点 数目 大 约 是 改进 


采样 点 数目 


184 个 、300 个 和 401 个 ; 与 三 次 Hermite 插值 拟 合算 法 
产生 的 欠 冲 点 数目 相 比 ， 


分 别 少 了 128 个 、249 个 、375 个 和 


与 经 


490; 


包 络 拟 合算 法 产生 的 欠 冲 点 数目 相 比 ， 分 别 少 了 


296 个 、555 


个 、860 个 和 1118 个 。 


显然 ， 本 文 算法 对 不 同情 绪 标 签 的 脑 电信 号 进行 
的 效果 优 于 经 


(database of emotion analysis using physiological signals, DEAP ) 


中 四 组 不 同情 绪 (两 种 正 性 情 
进行 了 包 络 拟 合 ， 计 算 经 


绪 和 两 种 负 性 情 


绪 ) 的 脑 
包 络 拟 合算 法 、 分 段 抛物 线 插 值 拟 


es a 
已 | 百 驴 


合算 法 、 三 次 Heimite 插值 拟 合 


算法 和 本 文 算法 在 包 络 拟 合 时 


产生 的 欠 冲 点 数目 。 结 果 如 表 1 所 示 。 

表 1 采样 点 为 2000 时 不 同情 绪 标 签 的 欠 冲 点 数 

算法 love exciting depressing hate 
经 典 拟 合 100 89 98 93 
分 段 抛物 线 64 55 46 47 
三 次 Hermite 62 41 45 32 
本 文 算法 30 19 17 18 

表 2 采样 点 为 4000 时 不 同情 绪 标 签 的 欠 冲 点 数 

算法 love exciting depressing hate 
经 典 拟 合 180 169 193 185 
分 段 抛物 线 121 110 102 88 
三 次 Hermite 107 90 85 74 
本 文 算法 60 41 36 35 

表 3 采样 点 为 6000 时 不 同情 绪 标 签 的 欠 冲 点 数 

算法 love exciting depressing hate 
经 典 拟 合 268 271 273 296 
分 段 抛物 线 161 165 146 151 
三 次 Hermite 143 145 128 132 
本 文 算法 84 52 46 66 


包 络 拟 合 
分 段 抛物 线 插值 拟 合算 法 和 


的 包 络 拟 合算 法 、 


三 次 Hermite 插值 拟 合算 法 。 


接着 使 / 


改进 的 包 络 拟 合算 法 进行 EMD 分 解 ， 以 分 解 的 


前 两 个 IMF 分 量 为 例 ， 
用 改进 的 包 络 拟 合算 法 进行 EMD 分 解 ， 每 一 个 IMF 分 量 能 
很 好 地 贴近 原始 信号 。 由 以 上 仿真 结果 可 知 ， 
的 包 络 拟 合算 法 ， 显 著 的 减少 了 欠 冲 点 数目 。 而 且 


结果 如 图 4 所 示 。 从 图 中 可 以 看 出 ,使 


本 文 提出 的 改进 
于 增 大 了 


极 值 序 列 ， 使 得 包 络 拟 合 的 曲线 ， 更 加 贴近 原始 信号 。 


xlt 


IMF1 
[=] 
pe 


100 120 140 160 180 200 


140 160 180 200 


图 4 EMD 分 解 的 结果 


4 ”结束 语 


本 文 针对 非 平 稳 随 机 的 脑 电 信号 ， 采 
的 方法 , 增加 了 极 值 点 数 


包 络 拟 合 产 4 


仿真 结果 证 明了 本 文 改进 的 包 络 拟 合算 法 的 适 


En 


用 扩大 极 值 点 的 范 
， 有效 地 减少 了 EMD 分 解 过 程 中 ， 

E 的 欠 冲 点 数目 ,得 到 的 包 络 线 有 着 更 好 的 平滑 性 。 
用 性 和 有 效 性 。 


但 本 文 算法 ， 


a) 经 验 模 态 分 解 的 “端点 效应 ”021。 


还 有 以 下 两 个 方面 的 问题 需要 解决 : 


由 于 信号 在 站 


点 处 在 


大 多 数 情况 下 不 是 极 值 点 ， 因 此 经 验 模 态 分 解 在 进行 包 络 线 的 


插值 拟 合 时 ， 


在 信号 的 两 端 会 产生 发 散 现 象 。 


同时 随 着 经 验 模 


态 分 解 的 不 断 筛 选 ， 引 起 较 大 的 失真 。 


录用 稿 


b) 求解 “ 伪 极 值 点 ”的 效率 。 在 对 “ 伪 极 值 点 ”进行 判断 
时 ， 对 信号 进行 了 裔 历 ， 因 此 效率 不 高 。 在 后 续 研 究 工 作 中 ， 
将 先 找 出 欠 冲 点 所 在 的 区 间 , 在 欠 冲 点 区 间 上 进行 “ 伪 极 值 点 ” 
的 判断 。 
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